Inverse slow relaxation in granular hopping systems 
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We present experimental results that demonstrate a glassy behavior in the conductance of quench 
condensed insulating granular metals that is different from that observed in continuous disordered 
systems. Exciting the granular system by biasing the sample with a high electric field results in a 
slow conductance change both during the excitation and during the relaxation back to its steady 
state. The time scales for these processes are many orders of magnitudes longer than the typical 
hopping time. We find that, initially, this conductance change has an opposite sign to that observed 
in similar experiments performed on continuous films. Only after relatively long times this trend 
is reversed and the samples exhibit conventional behavior. We suggest that the granular systems 
exhibit a new glassy process related to charge redistribution among the grains. This process combines 
with the relaxation processes that are characteristic of disordered systems to give rise to a unique 
relaxation profile. 

PACS numbers: 73.23.-b; 73.40.Rw; 75.50.Cc 



Glassy behavior of the conductivity of strongly disor- 
dered systems has attracted attention for over a decade. 
Exciting a hopping system out of its thermal equilibrium 
leads to an increase in conductivity. The relaxation pro- 
cesses toward equilibrium are characterized by extremely 
long times, memory phenomena and aging effects. The 
majority of experiments were performed on Indium Ox- 
ide thin films which were excited by applying a gate 
voltage, thus abruptly changing the carrier concentration 
[3, UTS 0, HI- Following the sharp initial increase, the 
conductivity exhibited a logarithmic time decay to a new 
equilibrium state with time scales of the order of many 
hours or days. This slow conductivity relaxation showed 
many phenomena characteristic of glasses such as scaling 
behavior in aging experiments and memory effects. 

Very similar glassy behavior (e.g., slow relaxation, 
memory and aging) was reported also in granular Au 
or Al [y, 0] demonstrating that this phenomena is not 
unique to Anderson insulators. This is not surprising 
since, in principle, the electrical transport in a granular 
metal in the insulating phase is very similar to that of a 
disordered continuous metal. The main difference is that 
in a granular system the electrons are localized on iso- 
lated grains instead of localized electronic states created 
by the random potential. Each grain is a small metal and 
can contain a large number of electrons, unlike the single 
localized state. 

Recently, Orlyanchik and Ovadyahu [H used a different 
technique to excite Indium oxide samples in the hopping 
regime. They biased the sample by applying a large elec- 
tric field, thus driving the system out of its ohmic regime, 
and studied the time-dependent conductance both dur- 
ing excitation (while applying the voltage) and relaxation 
(after the field was removed). They found that the ap- 
plication of the field caused a fast conductance increase 
followed by a slow rise towards an asymptotic value. Re- 
moving the field caused a fast conductance decrease fol- 
lowed by a very slow drop toward the thermal equilib- 



In this letter we report on an experiment in which 
we apply a similar technique to that of Orlyanchik and 
Ovadyahu to excite a granular metal. Surprisingly (un- 
like the case of applying a gate voltage), we find that 
applying a bias voltage gives rise to an opposite trend 
to that seen in the continuous disordered systems. Wc 
attribute this difference to the unique electric field dis- 
tribution in granular metals. 

The samples studied in this work were various (Ag, Au 
and Ni) thin granular films fabricated using the quench 
condensation technique. This method enables growing 
and measuring a 2D granular sample at low temperature 
and high vacuum without exposing it to the atmosphere 
or thermally cycling it [^, 1(J 11 1 . Prior to the cool-down, 
two Gold leads were evaporated on a insulating Si/SiO or 
glass substrate. The substrate was then placed in a UHV 
system and immersed in liquid He. An ultrathin metal 
film was then evaporated on the cold substrate, form- 
ing a discontinuous, granular structure with grain sizes 
of 100-200 A and intergr ain distances of a few tens of A. 
Sample resistance was monitored during the growth and 
the process was stopped at a desired resistance. We stud- 
ied samples with resistances in the range of 50 MCI — 5 
Gtt. After terminating the film deposition the sample was 
allowed to relax for a few hours until the conductance 
change was smaller than 0.1 %/hour and a base conduc- 
tance line was determined. Then, a large electric field, 
F n = 2 — 600 V/cm, was applied abruptly and the time- 
dependent-conductance, G(t), was monitored. After a 
desired waiting time, t w , the field was abruptly reduced 
to Fo k 0.5 V/cm and the low voltage relaxation G(t) 
was measured. All conductance measurements were per- 
formed by measuring the current through a small series 
resistor (at least two orders of magnitude smaller than 
the sample resistance), similar to the method described 
in ref [|. 

Typical excitation and relaxation G(t) profiles for two 
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FIG. 1: Relative conductance change versus time for two Ag 
granular films excited with a bias field ^=200 V/cm. (a) 
relaxation curve under F =0.5 V/cm, R=360 Mf2, T=4.2 K. 
(b) the relative conductance during the excitation, R=80 Mtt. 



different samples are shown in fig. 1. It is seen that the 
initial trend was opposite to that observed by Orlyanchik 
and Ovadyahu During the excitation the conductiv- 
ity decreased with time instead of the increase observed 
in continuous films, and during the relaxation process 
the conductivity increased slowly contrary to decrease 
observed in continuous films. After a period of time that 
depended on t w (typically, a few hundred to a few thou- 
sand seconds), the trend reversed and the conventional 
behavior typical of Anderson insulators was observed. 
This behavior was observed in all investigated granular 
metals without significant dependence on the material. 

Since the behavior during excitation and relaxation 
is similar (with opposite signs) for the reminder of this 
paper we will focus on the processes of relaxation af- 
ter switching off the electric field. In figure 1 we plot 
the change in conductivity versus time, AG(t)/G = 
{G(t) — G )/Go, where Go is the extremum conductance 
at which the change in the trend occurred and t=0 when 
the excitation was introduced or removed. We define two 
parameters, (as seen in fig. la): r, the time it takes the 
conductivity to reach G = Gq and A, the change in con- 
ductivity during r. In all measurements we identify two 
regimes of the relaxation G(t) profiles: the initial con- 
ductance increase at t < t we name the "opposite effect" 
and the conductance decrease at t > t which is the con- 
ventional behavior similar to the results in continuous in- 
sulators. We note that a typical sample had a resistance 



FIG. 2: (a) Conductance during relaxation for different val- 
ues of t w (b) The same data as in (a) plotted as a function 
of t/t w . R=5.4 GQ,, F„ = 100 V/cm. The conductance os- 
cillations which appear in the graphs are rather reproducible 
and probably result from the mesoscopic nature of the sample 



of ~ 100 Mfl and the measured circuit capacitance was 
~ 100 pF yielding RC ^10 ms, which is many orders of 
magnitude smaller than the typical measured relaxation 
time t. 

In order to study the glassy nature of the " opposite ef- 
fect" we performed aging experiments using the following 
scheme: We stabilized the system at low applied electric 
field, Fo, for several hours, thus reaching a steady state 
of the system. Then we excited the system by applying 
high bias field, F n , for a waiting time, t w after which 
the field was reduced back to F . Conductivity measure- 
ments were performed for a time of the order of 3 x t w 
after removing the high field. The relaxation curves of 
the conductivity versus time for three different t w are 
plotted in fig. 2(a). The same data is plotted in fig. 2(b) 
as a function of t/t w . It is seen that all curves collapse 
onto a single curve as a function of t/t w giving rise to 
simple scaling behavior [i. 13|. 

Varying the strength of the exciting field also had a 
significant effect on the G(t) profiles. Fig. 3a which 
depicts AG/G for different F n s shows that the relative 
conductance change, A, grows linearly with F n . For high 
enough F n , switching off the field gave rise to a short 
period of " negative conductance" (see fig. 3b), indicating 
that, for short times, the current flow was opposite to 



3 



o.o 



(9-2.0x10 

(9 

< 



-4.0x10 
5.0x10 s 



0.0 




(a) 



1 0OV/cm 
250V/cm 
400V/cm 
600V/cm 



F (v/cm) 



50 100 150 200 250 300 350 400 



C9 



-1.0x10 



-1.5x10 



(b) 



8000 



12000 16000 20000 

t(s) 

FIG. 3: (a) Relaxation curves after excitations by different 
F„. i m =2u00 sec and R=150 Mfi. The inset shows the am- 
plitude, A, as a function of applied bias voltage for a 320 Mil 
sample . (b) Relaxation processes where "negative conduc- 
tance" was obtained, R=250 MQ. 
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FIG. 4: The excess conductance versus time after switching 
off a 2 Hz AC field of 25 V/cm. 



the voltage direction. We found that r does not depend 
on the value of the exciting field and no correlation was 
found between A and t w . 

All the above results were obtained employing a DC 
voltage excitation. We performed similar measurements 
using an AC voltage excitation and measured the AC 
current through the sample. Fig. 4 shows that exciting 
the samples with an AC voltage, even with frequencies 
as low as 2 Hz, caused the "opposite effect" to disappear, 
leaving only the conventional, Orlyanchik and Ovadyahu 
relaxation trend. 

In the following paragraphs we propose a model for 
the observed results, based on the granular nature of the 
samples. The conductivity in these systems is governed 



by electron hopping between the metal grains 14[ . When 
an external electric field F n is applied to the system the 
charge on each grain reorganizes itself in order to reach 
electrostatic equilibrium. If the metallic grains are far 
from each other the charge on the grains is fixed and 
is redistributed on the grain so that the internal field 
cancels out the external F n . This process occurs over a 
characteristic metallic screening time (typically less than 
1 ps). However, for closely packed grains, strong elec- 
trons interactions cause the rearrangement of the charge 
distribution in each grain to affect the arrangement of its 
neighbors. Therefore both the magnitude and the direc- 
tion of the screening internal field in each grain change in 
time, mostly by hopping events, until it reaches a steady 
state. This process is interactive - hierarchic, hence is 
characterized by very long times leading to a slow de- 
crease of the electric current in the system. After the 
excitation is removed, the grains discharge in order to 
adapt to the new electrostatic potential. This too is a 
slow process leading to a long-time current decrease, in 
which the direction of the current is opposite to that 
of the original excitation. This leads to a new glassy 
process which competes with the conventional relaxation 
processes of a hopping system relaxing to equilibrium. 

In order to test this model we simulated the time de- 
pendence conductivity in a one dimension granular sys- 
tem. We consider a chain of metallic grains connected to 
two large leads which are biased by a voltage, V. Each 
grain is connected to its neighbor by junctions of capaci- 
tance Ci and tunneling resistance Ri , as illustrated in the 
inset of fig 5. The potential, , of each grain is given by 



Cjifij + n a e 



(1) 



where j is the index of the neighboring capacitors and 
grains, rij is the number of the electrons added to the 
grain i from the initially neutral condition and cs is given 
by c j ■ The potentials of the two leads are and V. 
The total electrostatic energy of the system is given by: 



U 



(2) 



By subtracting the total coulomb energy of two dif- 
ferent combinations of nj and V (state 1 and state 2), 
we find AE = U<2 — Ui, the total coulomb energy cost 
of the system transferring from one state to the other. 
However, the charge distribution in the system changes 
in order to keep a constant bias voltage, therefore we add 
to AE the energy supplied by the voltage source when 
an electron tunnels through any capacitor. 

We define a vector ~jt(V, ni, n*, t) which is the prob- 
ability to have n\....ni electrons on grains l....i respec- 
tively, under bias voltage V at time t. 

At non zero temperature, there is a probability of tun- 
neling processes, for which the rate of change from one 
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FIG. 5: Simulation of the current versus time, normalized 
to the hopping time, after an excitation with bias voltage 



V=500 V of a 5 grains chain, d 



10" 



F, Ri = 10 s MO. 



The inserts show a schematic description of the model and 
the dependence of the characteristic relaxation time on the 
number of dots in the chain. 



distribution to another is: 



TV„1 12 2^ 

1 [n 1 ...n i , n 1 ...n i ) 



1 



i-. 



-AE 



:), (3) 



Re 2 " I -exp{AE/k B T)' 

where n\...n 1 l is the charge distribution in the state k and 
R is the resistance of the junction. 

The master equation for the system is thus given by: 



d~P{V,ni, ■■■■,n i ,t) 
dt 



= T--p(y,n 1 ,....,n i ,t), (4) 



where T is the rate matrix and ~J? is the state vec- 
tor. Because the system is not in equilibrium, the cur- 
rent through each junction is not the same. The mea- 
sured current is the sum of all the possibilities to tun- 
nel through the last junction. Our initial conditions are 
probability 1 for zero electrons on each grain. We then in- 
troduce a high bias voltage V for a time of t w . We solve 
eq. 4 and calculate the current for every step. Fig. 5 
shows the relaxation curve of the current after switching 
off the high voltage excitation so that V=0. The current 
is negative and has an amplitude which decreases to zero 
over a time which is much longer (~3 orders of magni- 
tude) than a typical hopping process even for a system of 
5 grains. The inset of fig 5 shows that the characteristic 
relaxation times increase monotonically as a function of 
the number of grains, hence the times can be expected 
to be much longer in a macroscopic sample. We note 
that this is clearly a simplistic model which takes into 
account only nearest-neighbor hops and doesn't include 
the disorder properties of the system but we believe that 
it conveys the essential physics of the slow, glassy dy- 
namics responsible for the time dependent conductivity. 

In the actual experiment the electric field is not re- 
duced to zero but rather to a small dc field for conduc- 
tivity measurements. So the negative current observed in 
the simulations has to be added to a finite positive cur- 
rent. The sum of those two can be expected to lead to a 



current profile similar to those seen in figure 3b in which 
the total conductivity may begin negative but eventually 
approaches a steady state value of G . The change in the 
conductance is obtained as a response to the DC applied 
field. When the direction of the field changes faster than 
the redistribution process (i.e. as a result of an AC bias) 
the opposite effect does not show up. 

In summary, we presented results of the time depen- 
dent conductivity of insulating granular metals after exci- 
tation with a high electric field. We find that a disordered 
granular system behaves differently than its continuous 
counterpart. We interpret the results as originating from 
a competition between two different trends. The first 
process is the relaxation of the systems toward equilib- 
rium and is typical of all hopping systems. The second 
is a consequence of charge rearrangement on the metal- 
lic grains and is expected only in granular systems. In 
a sense, the granular metal excited by an electric field 
is similar to spin glass subject to a sudden magnetic 
field change. Despite the fact that a spin flip occurs on 
times of a few ns, the magnetic interactions in the system 
cause the arrangement of the whole system to occur on 
timescales which are many orders of magnitude larger. 
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